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Abstract. We study the steady state behaviour of a confined quantum Brownian 
particle subjected to a space-dependent, rapidly oscillating time-periodic force. To 
leading order in the period of driving, the result of the oscillating force is an effective 
static potential which has a quantum dissipative contribution, Vqd, which adds on 
to the classical result. This is shown using a coherent state representation of bath 
oscillators. Vqd is evaluated exactly in the case of an Ohmic dissipation bath. It is 
strongest for intermediate values of the damping, where it can have pronounced effects. 
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1. Introduction 

Parametric phenomena, induced by giving a time-dependence to a parameter or 
coupling constant, are ubiquitous in physics. Applications range from confining ions 
in quadrupolar traps (Paul traps) [lj to controlling particle bunching and dilution in 
particle accelerators [2]. 

When the parameter in question oscillates rapidly, certain simplifications occur, as 
first demonstrated by Kapitza for the case of a simple pendulum whose point of oscil- 
lation is vibrated rapidly [3]. Kapitza's treatment for the pendulum was generalized 
by Landau and Lifshitz to apply to any forced system evolving through Newtonian dy- 
namics, provided the forcing depends on the spatial coordinate and oscillates rapidly 
enough in time [U [5]. The result is simple and elegant. To leading order in the pe- 
riod of driving, the system feels an additional effective static potential; the form of this 
Kapitza-Landau-Lifshitz (KLL) potential, Vkll, depends on the spatial variation of 
the forcing. 

Recently, this treatment has been generalized to the case of a Brownian particle 
subjected to similar forcing 0,[7]. The result of classical dissipation (CD), to leading 
order, is again an effective static potential, but with additional Wronskian terms, Vcd, 
reflecting the effect of dissipation in this system. To next order, it is shown that the 
result cannot be written as an effective potential, except in the limit of high damping 
[6]. An approximate Langevin approach used in [TJ, leads to a result which agrees with 
the correct one to leading order, and to the next order if the damping is very large. 

The purpose of the present paper is to address this problem for a quantum dissipa- 
tive system. In the absence of dissipation, Cook et al. [12] showed that to leading order 
a Schrodinger equation with a rapid periodic driving term is tantamount to one with 
an additional static potential of the KLL form, while Rahav et al used an expansion 
in powers of the period to show that to the next order, the effective static Hamiltonian 
involves the particle coordinate and momentum in a non-separable way [13]. In this 
paper, we model quantum dissipative (QD) Brownian motion by considering a particle 
coupled to an infinite set of harmonic oscillators [H [H [10l [11] . We are interested in 
steady state properties, which are independent of the initial condition. Using an initial 
coherent state of bath oscillators, a c-number generalized quantum Langevin equation 
(GQLE) is derived. We show that to leading order, the rapid forcing is tantamount to 
an effective static potential which can be written as Vkll + Vcd + Vqd- Here Vqd is 
an explicitly quantum dissipative contribution to the potential, which vanishes in the 
classical limit [6]. In the complete absense of dissipation, the result reduces to Vkll in 
agreement with [121 HB]- Vqd has most pronounced effects at intermediate values of the 
damping. In that regime, Vqd is found to strongly influence the shape of the potential, 
e.g. the number and nature of the extrema of the potential. 
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2. Model and Formalism 

We consider a particle coupled to a model heat bath of the Zwanzig form [8], described 
by the following Hamiltonian : 

H = f- + V {x) + V^{x,t) + E[t^- + J^jfe - -^M)"], (1) 
2m j l 2mj 2 J v 

where x and p are coordinate and momentum operators of the Brownian particle, {qj,Pj} 
is the set of co-ordinate and momentum operators for the heat-bath oscillators, and Vq{x) 
is the harmonic confining potential. The potential V ext (x, t) arises from the external force 
field. The co-ordinate and momentum operators follow the usual commutation relations 

[x,p]=ih, [q j ,p k ]=ih5 jk . (2) 

Eliminating the reservoir degrees of freedom [9j [11] , we obtain the operator GQLE 
corresponding to Hamiltonian (1) : 

mx(t) + m ( dt xl {t - ti)x(tx) + V£(x) = F ext {x, t) + f)(t), (3) 
J o 

where noise operator fj(t), memory kernel j(t) and external force F ext (x,t) are given by 
m = £ c 4fe(°) - -^-a^CO)) cos(^t) + M2i sin(a; t )l (4 ) 



7W=E- 2 3COs(^), (5) 

UM) = -C(M). (6) 

Our aim is to replace the exact GQLE in operator form Eq. (3) by its c-number 
equivalent. To achieve this, we follow the method of Banerjee et al. [TH [To] and 
introduce product separable quantum states to describe the initial t — state of the 
particle and the bath oscillators, 

\lpinMal) = |0){K-}}, (7) 

where \<f>) denotes any arbitrary initial state of the particle and |aj) corresponds to the 
initial coherent state of the jth oscillator. Here \aij) is given by 



ai> = «p(-^f) E(-rr)l^>. ( 8 ) 



n,=0 



and ctj is expressed in terms of the variables of the coordinate and momentum of the 
jth oscillator 



(&(0)) - <x(0)> = \/^(«* + «i) (9) 



<PM = iy^r(<*s-<*i)- ( 10 ) 



Quantum Brownian motion under rapid periodic forcing 



4 



Following [T4J [T2] and doing the quantum-statistical averaging starting from the intial 
product separable quantum state (Eq. 7), one obtains the c-number GQLE : 

mx(t) +m f dti7(t - t x )x{t x ) + Vq(x) + V^ xt (x, t) = n(t) + Q (x, t) 
Jo 

+ Qext{x,t), (11) 

where 

x{t) = (x), 



V(t) = (V(t)) = E[{(?i(°)) - (x(0))}-^cos(^) + /^(^.(0))sin(^), (12) 

j m j UJ j V 3 j 

and Qo and Q ext are quantum fluctuation terms, given by : 

Qo(x 1 t) = VZ(x)-(VZ(x)) 1 (13) 

Q ext (x,t) = V^{x,t) - (K* t 0M)>- (14) 
In order that rj(t) be an effective c-number noise, we must have 

(V{t))s = 0, (15) 

(V(t)v(t')) s = iE^4^ coth (^)cos^(t-t'), (16) 

J J J 

where (....) s denotes statistical average over the initial distribution of the mean values 
of the momenta and co-ordinates of the bath oscillators. Equations (15) and (16) 
imply that £ (t) is centered around zero and satisfies the quantum fluctuation-dissipation 
relation, and are obtained if and only if the initial quantum mechanical mean values of 
momenta and co-ordinates of the bath oscillators have the following distribution [141 [T5] : 

where fij = [exp(|^.) — l] -1 , the average thermal phonon number at temperature 
T. Thus, the statistical average of any quantum mechanical mean value 
Oj((Pj{0)), {(qj(0)) - (x(0))}) is defined as 



(Oj), = J OH<Pi(0)>, {<^(0)> - (m)})Pjd<Pj{0))d{W)) ~ (I 8 ) 

Using equations (12), (17), and (18), one can easily show the properties (15) and (16) 
of the c-number noise. Pj is a canonical Wigner distribution for a displaced harmonic 
oscillator and always remains positive [16]. Now, one can interpret equation (11) as a 
c-number GQLE which is governed by a c-number noise r)(t) originating from the heat 
bath characterized by the properties (15) and (16). The two quantum fluctuation terms 
Qo(x,t) and Q e xt(%,t) originate from the nonlinearity of the confining potential V (x) 
and the externally applied rapidly oscillating potential V ext (x, t), respectively. 
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In this section, we discuss the derivation of the quantum terms Qo(X, t) and Q ext (X, t). 
Based on the quantum nature of the system, one can write 

5x(t)=x(t)-x(t), (19) 
Sp(t)=p(t)-p(t), (20) 

where x(t) = (x(t)), pit) = (pit)). Then, making a Taylor series expansion around x in 
Eqs. (13) and (14), one obtains 

Q (x,t) = -EV +1) W(^">. (21) 

n>2 U - 

Q ex t(x,t) = - £ -V}: t +1 \x,t)(8x n ), (22) 

where V {n+1 \x) and v£ +1) (a, t) are the (n+l)th derivatives of Vq(x) and V ext (x, t). The 
terms Qo(x, t) and Q e xt(x,t) depend on the quantum factors (5x n (t)) whose evaluation 
is discussed below. 

Substituting Eqs. (19) and (20) into Eq. (3), we obtain 

m5x + m f dta(t - h)5x + Vj(x)6x + V? xt (x, t)5x + ]T ^-V^ n+1 \x)(5x n ) 

-* n>2 U - 

+ E -y e ( : t +1 \x,t)(5x n ) = 5m, (23) 

n>2 n - 

where 5fj(t) = rj(t) — r)(t). 

We now consider the case of a harmonic confining potential V (x) = ^uiuJqX 2 , in 
which case the derivatives in Eq. (21) vanish, implying Qo(x, t) = 0. With rapid driving, 
excursions around the smooth part are small, so that for our purpose of deriving the 
leading order term, we may replace (5x 2 (t)) by {Sx 2 (t)) where (....)o denotes an average 
in the absence of driving. Thus, we have 

Q ext (x,t) ~ -ly e -(x,0(^ 2 (0)o. (24) 

Equation (23) then reduces to 
..rt 

mbx + m / dti'jit — ti)5x + 17110^5 x = S^. (25) 
Jo 

Equation (25) may be solved by Laplace transformations, leading to 

6x(t) = G(t)6x(0) + H(t)5x(0) + /* dt'Hit - t'Wit'), (26) 

Jo 

where G(t) and H(t) are the inverse Laplace transforms of G(s) and H(s) respectively 
given by 

^, s s + ^(s) 
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Here 



-f(t)e- st dt 



(29) 



is the Laplace transform of the Frictional kernel j(t). Equation (26) then leads to 



(Sx 2 (t)) = G 2 (t)(Sx 2 (0)) + 



H\t) 



(^ 2 (0)) + 



G(t)H(t) 



ni 



(5x{0)5p(0) + 5p(0)Sx(0)) 



+ 2 I* dt' f dt"H(t - t')H(t - f)(6fj(tf)Sfj(f)). 
Jo Jo 



(30) 



Let us choose initial conditions corresponding to minimum uncertainty states 
[HI9]:, so that (6x 2 (0)} = (^ 2 (0)> = and {6x(0)Sp(0)+Sp(0)6x(0))o = K. 

In order to find H{t) and Git) (or equivalently H(s) and G(s)), we need to know ^(t). 
We make the customary choice of an Ohmic heat bath, which leads to 7(i) = jo5(t) and 
finally to an explicit forms for H(t) and G(t). In the underdamped regime (uq > 7/2), 
we find 



Hit) = —e~ 
G(t) = e"^ 



2 r sm(o;it), 
cos(a>i£) + 



7o 



sin(a>ii) 



(31) 
(32) 



where iO\ = ±y uj 2 — ^f. For the overdamped case (lu < io\ becomes imaginary and 
H{t) and Git) are modified to 



H(t) = — e^'smli^t), 



G(t) 



_2a+ 
e 2 



coshfa^i) H -^-r sinhfu/i) 

2muj l 



(33) 
(34) 



where u[ = ± y ^ — co>o- 

With the Ohmic condition for the heat bath, 7(i) = 7 <5(t), the double integral in 
equation (30) can be evaluated with the result 

2 



7o 
m7r Jo 



dujhu coth 



2k B T 



1 - e-^-^fcoswit + (f - ^i) 2 ^] 



uj 2 - + ^7o^ 



It reduces to a simple form for time large compared to 7 1 : 



(^ 2 )o 



To 



[ dujfiLU coth ( \ 

Jo ^2k B T> 



^-u 2 y + 1 W 



In the limit of weak damping (70 << Uo), Eq. (35) reduces to 
(Sx 2 ) 



-A_ coth 

while with strong damping (70 >> cjq), we obtain 



(^ 2 >c 



h 



7T70 



21n( 



7o 



, . mhuo 2 
+ H: 



mhuj 2 ^ 



+ IE 



(35) 



(36) 



(37) 



~mu ' '27T7 /cbT mhuj 2 ^ h^ 
where 7^ ~ 0.57772 is the Euler's constant. When substituted in Eq. (24), these 
answers for {Sx 2 )q determine the leading order form of the quantum fluctuation term 
Eq. (22). 
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4. The effective potential 

In this section we show that on time scales larger than the period of driving, the 
leading order efffect is to produce an effective potential which has an additive quantum 
dissipative term Vqd whose form is derived below (Eq. 50). 

For the harmonic confining potential and Ohmic dissipative bath, we have seen that 
the operator GQLE reduces to 

mx(t) + rwy x(t) + Vq(x) = F ext (x, t) + Q ext {x, t) +r](t). (38) 

If r](t) is delta correlated as it is for an Ohmic dissipative bath, one can show that the 
noise term does not enter the effective Hamiltonian to leading order [20] . Now, following 
the KLL approach pE] one can write 

x(t)=X(t)+aX,t), (39) 

where X(t) and £(X, t) are the slow coordinate and rapidly oscillating coordinate 
respectively. Inserting Eq. (39) in Eq. (38), we obtain 



mX{t) + mf = - 7 X - 7o£ 



dV (X) d 2 V 



dx * ox 2 

+ F ext (X, t) + tF' ext {X, t) + Q ext (X, t) + & ext (X, t). (40) 
Equation (40) contains both slow and fast part and they must be separately equal. Thus 
m£ + l0 £ = F ext {X,t). (41) 
For F ext (X, t) = f(X) cos(fit) + g(X) sin(fit), we may choose 



1 



m{Q 2 



(f(X) + ^-g(X))cos(nt) 
mil 



7o 



+ (g(X)--^f(X))sm(Qt) 



(42) 



Now, putting Eq. (42) in Eq. (40) and averaging over one time-period (which is denoted 
by subscript r) of the external field, we obtain 

dVo(X) 



mX(t) + 7o X 



dX 



+ m xt (X, t)) T + (Q ext (X, t)) T + (ZQ' ext (X, t)) T . (43) 



Now, 



m xt (x,t)) T 



4m{Q 2 + 



^(f(X)+g 2 (X)) 



27o 



(Qext(X, t))i 



+ 



2T Jo 



(f'(X)g(X)-g'(X)f(X)) 



dt{f"\X) cos(ftt) + g"\X) sin(fit)} (5x 2 ) i 



(44) 



(45) 



Since the term in square bracket is independent of time (see Eq. (35)), we obtain 
(Q ex t(X, t)) T = 0. In the same manner one can show that in the limit of large t > 'Jq 1 

(£Q' ext (X,t)) T = 1 



+ 



2m(W + ^ 
^(f"(X)g(X)-g'"(X)f(X))} 



(5x 2 )} o {(f(X)f"(X)+g(X)g'"(X)) 

(46) 
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where we have replaced {5x 2 (t)) by its equilibrium value (5x 2 )q. Finally, one can write 
down the effective equation of motion for the slow variable as follows 

ml + 7oI = "^' (47) 



where 



V 



U eff = V + V cl + V QD , (48) 
1 r, o,„ss 270 



2 1 



■m- 



2 



(/ 2 (X) + 5 2 (X)) + ^ ^ dy(g{y)f'(3,) - /(y)<fly)) 



= V^ll + ^cd, (49) 

with V KLL = (/2(X)+g y )} and V CD = f x dy(g(y)f(y) - f(y)g'(y)). Finally, 

the quantum contribution is given by : 

' \(5x 2 )}\ [ X dy(f(y)f"'(y)+g(y)g"'(y)) 



Vr 



QD 



2 1 % 



in- 



2 , 



+ ^ [ dy(g(y)f"'(y) - f(y)g"'(y))}. (50) 

Equation (50) embodies the leading order quantum correction. In the classical limit, 
Vqd vanishes and the effective potential reduces to V + Vrll + Vcd, which is the result 
of [6j [7]. We also note that in the absence of dissipation 70 — > 0, both Vcd an d Vqd 
vanish and the effective potential reduces to Vq + Vkll, consistent with the results of 
[T2l [13] . Finally, in the limit of very high dissipation 70 — > 00, the quantum correction 
is negligible as is evident from Eq. (35). Thus, we conclude that the effects of quantum 
dissipative contribution are largest for intermediate values of damping. 

To illustrate the evolution of the form of V e ff as the damping is varied, we consider 
a system which is perturbed by an oscillatory potential with a Gaussian profile i.e. 
Vext{%,t) = rexp(— x 2 /a 2 ) cos(fit). The effective potential then has the following parts 

Vo ==^=X 2 , (51) 

fro 

V d = = ^ = ^-^X 2 exp(-2X 2 ), (52) 
-fro 1 + -Jl 



v QD = = ^ = - A 72 J (^ 2 )] [6(r(2, 2X 2 ) - r(2; 



+ 3(1 - e- 2 ^ 2 ) - (r(3, 2X 2 ) - T(3))] , (53) 

where T(b, y) is the incomplete gamma function and Eq = ^mu 2 a 2 . The above equa- 
tions are expressed in terms of dimensionless variables X = X/a, A = 16m ^ Q2 a 2 , and 
[(5x 2 )]o = [(Sx 2 )]o/a 2 . Figures (1), (2), and (3) depicts the effective potential profile in 
the low, intermediate and high damping cases respectively. As remarked in the previous 
paragraph, we see that the effects of Vqd are most pronounced at intermediate values 
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Figure 1. (color online) Plot of the effective time independent potential versus X for 
very weak damping with external driving potential V ext (x, t) —T exp(— x 2 /a 2 ) cos(fit) 
(a) without (red filled circle; Vo + V c i) and with (blue filled square; Vo + V c i + Vqd) 
quantum dissipative part, (b) only quantum dissipative part (Vqd) and (c) only 
classical part (V c i)- To plot this figure we used 70/mwo = 0.08, fl/ujo — 150, A = 17.85, 




Figure 2. (color online) Plot of the effective time independent potential versus X for 
moderate damping values with external driving and parameters as in Figure 1, except 
ja/mujQ = 1.0. 

of the damping. The presence of higher order derivatives of space dependent terms 
of forcing in the quantum dissipative part, Vqd, brings in some interesting changes in 
the qualitative form of the effective potential. It is evident from figure (1) that in the 
absence of dissipation the effective potential has two metastable minima and one stable 
minimum . But, on adding Vqd two metastable states become stable minima and the 
stable minimum becomes a maximum. 
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Figure 3. (color online) Plot of the effective time independent potential versus X 
for high damping values with external driving and parameters as in Figure 1, except 
ja/muiQ — 10.0. 



5. Conclusion 

In this paper, we have analyzed the properties of a quantum particle in a harmonic 
potential, coupled to a dissipative Ohmic heat bath, in the presence of a space depen- 
dent rapidly oscillating field. Starting from the usual heat-bath reservoir model, we 
obtained a c-number GQLE from the operator equivalent by using a coherent state 
representation of bath oscillators [HI US]- This enabled us to apply the Kapitza idea 
of separating into a slow part and a fast part, which consists of the rapid oscillations 
around the slow motion. The resulting equation for the slow motion is treated pertur- 
batively. To leading order, we find that the rapid forcing leads to an effective static 
potential which has a quantum contribution, Vqd, which adds to the classical result. 
This quantum contribution vanishes for both very small and very large values of the 
damping, but can have an appreciable and interesting effects for intermediate values. 
For instance, the number of stable minima in the effective potential can change, as il- 
lustrated by the case of a Gaussian oscillating force, for moderate values of the damping. 
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